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The  equilibrium  field  near  the  tip  of  a  crack  for  finite 
plane  strain  of  Incompressible  elastic  materials* 

by 

Rodney  A.  Stephenson 
California  Institute  of  Technology 


This  investigation  is  concerned  with  the  deformations  and  stresses 
in  a  slab  of  all-around  infinite  extent  containing  a  traction-free  plane 
crack,  under  conditions  of  plane  strain.  The  analysis  is  carried  out 
within  the  framework  of  the  fully  nonlinear  equilibrium  theory  of  homo¬ 
geneous  and  isotropic  incompressible  elastic  solids.  For  a  fairly  wide 
class  of  such  materials  and  general  loading  conditions  at  infinity,  as- 
symptotic  estimates  appropriate  to  the  various  field  quantities  near  the 
crack-tips  are  deduced.  For  a  subclass  of  the  materials  considered,  these 
results  —  in  contrast  to  the  analogous  predictions  of  the  linearized 
theory  —  lead  to  the  conclusion  that  the  crack  opens  up  in  the  neighborhood 
of  its  tips  even  if  the  applied  loading  is  antisymmetric  about  the  plane 
of  the  crack,  (e.g..  Mode  II  loading).  It  is  shown  further  that  the  non¬ 
linear  global  crack  problem  corresponding  to  such  a  loading  in  general 
cannot  admit  an  antisymmetric  solution. 

_ 

The  results  communicated  in  this  paper  were  obtained  in  the  course  of  an 
investigation  supported  in  part  by  Contract  N00014-75-C-0196  with  the 
Office  of  Naval  Research  in  Washington,  D.C. 


Introduction 


Early  investigations  of  crack  problems  beyond  the  scope  of  the  clas¬ 
sical  theory  of  elasticity  typically  retain  the  kinematic  assumption  of 
infinitesimal  deformations  while  relinquishing  the  linear  stress-strain 
law  in  favor  of  nonlinear  constitutive  relations  of  one  kind  or  another.^" 
Apparently  the  first  investigation  of  a  crack  problem  within  the  fully 
nonlinear  equilibrium  theory  of  elastic  materials  is  due  to  Wong  and  Shield 
[2].  In  [2]  an  approximative  global  solution  is  obtained  to  the  problem 
of  a  finite  crack  in  an  all-around  infinite  thin  incompressible  sheet  of 
a  neo-Hookean  material,  subjected  to  biaxial  tension  at  infinity,  on  the 
assumption  that  the  deformations  are  large  throughout  the  sheet. 

Knowles  and  Sternberg  [3],  [A]  deal  with  the  problem  of  an  infinite 
slab  of  a  compressible  elastic  material  containing  a  plane  traction-free 
crack,  within  the  equilibrium  theory  of  finite  plane  strain;  they  restrict 
their  attention  to  a  loading  of  uniaxial  tension  at  right  angles  to  the 
plane  of  the  crack  (Mode  I  loading)  and  analyze  the  structure  of  the  en¬ 
suing  elastostatic  field  near  the  crack-tips.  The  Mode  III  crack  prob¬ 
lem,  in  which  the  loading  at  infinity  is  one  of  longitudinal  shear  parallel 
to  the  edges  of  the  crack,  is  explored  asymptotically  in  the  finite  theory 
for  a  class  of  incompressible  elastic  materials  by  Knowles  [5].  Additional 
work  on  the  nonlinear  Mode  III  problem  may  be  found  in  [6],  [7],  [8].  The 
investigations  referred  to  above,  as  well  as  some  related  studies  in  finite 
elastostatics,  are  reviewed  in  several  survey  papers  [9],  [10],  [11],  [12]. 

^See,  for  example,  a  comprehensive  survey  by  Rice  [1]. 
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The  present  study,  which  is  also  set  within  finite  elastostatlcs, 
retains  the  body  geometry  underlying  [3],  [4],  as  well  as  the  assumption 
of  plane  strain,  but  supposes  the  material  to  be  incompressible.  The 
particular  constitutive  assumptions  Introduced  here  allow  for  both  harden¬ 
ing  and  softening  behavior  in  simple  shear,  but  preclude  a  loss  of  ellip- 
ticity  of  the  governing  displacement  equations  of  equilibrium.^  In  ad¬ 
dition,  more  general  loading  conditions  at  infinity  are  considered,  so  as 
to  encompass  in  particular  the  Mode  II  loading  of  simple  shear  parallel 
to  the  crack  faces. 

Section  1  contains  a  review  of  some  prerequisites  from  the  finite 
theory  of  plane  strain  for  homogeneous,  isotropic  Incompressible  elastic 
solids  and  introduces  a  particular  class  of  materials  underlying  the  sub¬ 
sequent  analysis.  In  Section  2  we  formulate  the  global  crack  problem 
with  which  we  are  concerned  and  recall  certain  properties  of  the  solution 
to  the  analogous  crack  problem  in  the  linearized  theory.  We  then  show 
that,  at  least  for  the  Mooney-Rivlin  material,  the  global  nonlinear  Mode  II 
problem,  in  contrast  to  its  counterpart  in  the  linear  theory,  does  not 
admit  a  solution  antisymmetric  about  the  plane  of  the  crack.  Section  3 
and  Section  4  are  devoted  to  an  asymptotic  analysis  of  the  elastostatic 
field  near  the  crack-tips.  Finally  Section  5  contains  a  discussion  of 
the  main  conclusions  reached  concerning  the  ensuing  local  deformation  and 
stresses. 


potential  loss  of  ellipticity  in  the  context  of  the  Mode  III  prob¬ 
lem  is  provided  for  in  [7],  (8]. 


-3- 


1.  Preliminaries  from  the  finite  theory  of  plane  strain  for  incompres¬ 


sible  elastic  solids 


Throughout  this  paper  we  shall  be  concerned  with  plane  finite 
elastostatics  of  homogeneous  and  Isotropic  Incompressible  materials, 
in  the  absence  of  body  forces.  Consider  a  body  which  in  an  undeformed 
configuration  occupies  a  cylindrical  region  of  space.  Denote  by  II 
the  open  projection  of  this  region  onto  a  plane  that  is  perpendicular 
to  the  generators  of  the  lateral  boundary.  Introduce  a  two-dimensional 
rectangular  cartesian  coordinate  system  in  this  plane  and  thus  let 
(x^,  x2)  be  coordinates  associated  with  the  position  vector1  x  . 

A  plane  deformation  of  the  body  parallel  to  the  (x^,  x2) -plane 


is  characterized  by  the  mapping 


y  =  y(x)  *  x+u(x)  on  II 


Thus  y  is  the  position  vector  of  points  in  the  deformation  image 


(1.1) 


n  * y1 (H) ,  while  u,  is  the  displacement  field.  We  suppose  the  function 
y  to  be  twice  continuously  differentiable  on  II,  and  the  mapping  (1.1) 
to  be  invertible. 

Let  F  be  the  associated  deformation-gradient  field,  so  that 


F «  v£  on  n 


(1.2) 


Since  the  body  is  incompressible,  the  deformation  (1.1)  must  be  locally 

^ere  and  in  what  follows,  letters  underlined  by  a  tilde  designate 
vectors  and  second-order  tensors  in  two  dimensions. 
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volume-preserving  and  hence 


Further,  let 


J  =  det  F  *  1  on  II  . 


(1.3) 


I-trFTF  on  II  .  (1.4)1 2 

With  the  aid  of  (1.3),  the  scalar  deformation  invariant  I  is  readily 
found  to  obey 

U2  on  n  ,  (1.5) 

and  1-2  only  in  an  undeformed  configuration. 

Next,  let  £  stand  for  the  two-dimensional  Cauchy  (true)  stress- 
tensor  field,  regarded  as  a  function  of  position  on  'II  .  For  an  equi¬ 
librium  deformation,  in  the  absence  of  body  forces, 

diVT  =  0,  t  =  on  JI*  .  (1.6) 

/v  ew  /v 

If  <j  designates  the  associated  Piola  (nominal)  stress  field  on  II,  one 
has  in  view  of  (1.3), 

o(x)  -  t(y(x))F  T(x)  for  all  x  in  IT  (1.7)^ 

In  terms  of  the  nominal  stresses,  the  equilibrium  equations  (1.6)  become 

divo-0,  a  FT  -  F  on  IT  (1.8) 


1The  superscript  T  indicates  transposition. 

2  -T  T 

Here  F  denotes  the  inverse  of  F  . 
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For  future  reference  we  also  recall  the  following  result.  Suppose 
T  is  a  regular  arc  in  II,  and  r**y(r)  its  image  in  II*.  Next  let  n 

rw  ^ 

and  n*  be  the  oriented  unit  normal  vectors  of  T  and  T*,  respectively, 
and  denote  by  s^  and  £  the  Piola  and  Cauchy  traction  vectors  defined  by 

s  *  a  n  on  T  ,  t  *  t  n*  on  T*  .  (1.9) 

<v  rw  /v  ' 

Then, 

£*£  on  T  if  and  only  if  £=  0  on  T*  .  (1.10) 

Moreover,  if  the  deformation  and  the  nominal  stress  field  are  suitably 
regular  in  the  closure  of  II,  (1.10)  continues  to  hold  true  for  an  arc 
T  on  the  boundary  of  II.  The  foregoing  result  thus  enables  one  to  im¬ 
pose  boundary  conditions  at  a  traction-free  edge  without  involving  the 
unknown  deformation  image  of  such  a  part  of  the  boundary. 

Suppose  now  the  elastic  solid  under  consideration  possesses  an 
elastic  potential,  and  let  W  stand  for  the  strain-energy  density  per 
unit  undeformed  volume,  regarded  as  a  function  of  position  on  II.  For 
plane  deformations  of  a  homogeneous  and  isotropic  incompressible  elastic 
solid,  W  depends  on  the  material  position  vector  exclusively  through 
the  deformation  invariant  I.  Hence, 

W(x)  ■  W(I(x) )  for  all  x  in  II  ,  (1.11) 

rw  rw  rv 

© 

where  W(I)  is  the  plane-strain  elastic  potential.  The  appropriate 
constitutive  law,  as  far  as  the  in-plane  stress  response  is  concerned, 


now  becomes 
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0*  2W' (I)F- pF-T  on  n  , 

<v  rv  ‘  /v  ' 


(1.12) 


in  which  W  is  the  derivative  of  W,  while  p  is  an  arbitrary  hydro¬ 
static  pressure  field  arising  from  the  constraint  of  incompressibility. 
In  view  of  (1.12),  the  second  of  the  equilibrium  equations  (1.8)  is 
satisfied  &  priori . 

We  now  turn  to  certain  restrictions  to  which  the  elastic  potential 
W(I)  is  subject.  First,  we  shall  take  for  granted  that 


W(I)  iO  for  I  >2,  W(2)  =  0  , 


(1.13) 


so  that  in  particular  the  strain-energy  density  vanishes  in  the  undeformed 


state.  Next,  according  to  the  relevant  Baker-Ericksen  inequality. 


W'  (I)  >  0  for  I  >  2  . 


(1.14) 


This  inequality  admits  a  simple  physical  interpretation  in  the  context 
of  a  plane  homogeneous  deformation  corresponding  to  simple  shear.  In 
this  case,  the  components  of  the  deformation-gradient  tensor  are  given  by 


1  k 


0  1 


(1.15) J 


in  which  k  is  the  amount  of  shear.  It  follows  from  (1.12),  (1.7)  that 
for  such  a  deformation  the  true  shearing-stress  component  obeys 


T12  =  t(k)  -  2kW'  (2  +  kZ)  , 


(1.16) 


LGreek  subscripts  have  the  range  (1,2)  and  summation  over  repeated  sub¬ 
scripts  is  implied. 
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and  we  henceforth  refer  to  the  graph  of  x(k)  as  the  response  curve 

for  simple  shear.  The  inequality  (1.14)  thus  holds  provided  x(k) 

and  k  have  the  same  sign,  as  demanded  by  obvious  physical  considerations. 

Finally,  we  recall  from  Abeyaratne  [13]  that  the  pertinent  displace¬ 
ment  equations  of  equilibrium  —  obtained  by  elimination  of  o_  between 
(1.8),  (1.12)  and  subject  to  the  constraint  (1.3)  —  are  elliptic  at  a 
solution  u  and  at  a  material  point  x  if  and  only  if 

W7  (1)  +  0  ,  2f,  ^ (I  -  2)  +  1  >  0  ,  (1.17) 

w'(I) 

where  I»I(x^  is  the  associated  value  of  the  deformation  invariant  de¬ 
fined  in  (1.4).  The  first  of  these  inequalities  is  assured  by  (1.14), 
while  the  second  is  readily  seen  to  be  equivalent  to  the  requirement 
that  the  slope  of  the  response  curve  for  simple  shear  be  non-negative  at 
k-  /I  -  2  . 

At  this  point  we  introduce  a  particular  class  of  incompressible 
elastic  materials,  hereafter  referred  to  as  power-law  materials,  which 
will  play  an  important  part  in  the  analysis  to  follow.  This  class  of 
materials  is  characterized  by 

W(I)  =  AIn+ BIn-1+ o(In_1)  as  I- oo  ,  (1.18) 

in  which  A,  B  and  n  are  constants  with  A  and  n  positive.  We 
shall  assume  further  that  the  asymptotic  identities  obtained  by  two  suc¬ 
cessive  formal  differentiations  of  (1.18)  are  also  valid.  It  should  be 
emphasized  that  (1.18)  merely  stipulates  the  asymptotic  behavior  of  the 
elastic  potential  at  large  values  of  I,  and  does  not  otherwise  restrict 
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its  specific  form.  A  special  case  of  an  incompressible  material  con¬ 
forming  to  (1.18)  is  supplied  by  the  Mooney-Rivlin  material,  with  the 
complete  plane-strain  elastic  potential 

W(I)«-|(I-2)  .  (1.19) 

This  particular  potential  evidently  satisfies  (1.18)  with  n=l,  A=p/2, 

B  =  -p . 

According  to  (1.16),  the  response  of  a  power-law  material  in  simple 
shear  obeys 


t (k)  =*  2nAk^n  ^+o(k^n  ^)  as  k-*oo  .  (1.20) 

The  corresponding  response  curve,  for  a  range  of  values  of  the  parameter 
n,  is  indicated  in  Figure  1.  For  n>l  the  material  hardens  under 
simple  shear,  for  n<  1  it  softens,  while  for  n=l  the  response  function 
t(k)  is  asymptotically  linear.  For  this  reason  we  shall  from  here  on  refer 
to  n  as  the  "hardening  parameter."  When  n=*l/2,  x(k)-*A  as  k-*oo; 

when  n<l/2,  x(k)  is  bound  to  reach  a  maximum  beyond  which  it  monoton- 
ically  declines  to  zero.  Such  a  "collapse  under  shear"  is,  as  previously 
remarked,  associated  with  a  possible  loss  of  ellipticity  of  the  governing 
field  equations.  In  contrast,  if  n>l/2,  the  ellipticity  conditions  are 
necessarily  met  at  all  sufficiently  large  values  of  I. 
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2.  Formulation  of  the  global  crack  problem.  Global  considerations. 

We  proceed  here  to  the  global  formulation  of  the  problem  with  which 
we  are  concerned.  Let  £  be  the  straight-line  segment 

£={x|x2!=0,  -bix^b]  ,  (2.1) 

and  suppose  II  to  be  the  region  exterior  to  <£,  as  indicated  in  Figure  2. 
Accordingly,  the  undeformed  body  is  an  all-around  infinite  slab  containing 
a  plane  crack  of  length  2b.  We  shall  assume  that  the  faces  of  the  crack 
are  traction-free  and  that  a  known  homogeneous  state  of  plane  deformation 
prevails  at  infinity. 

The  problem  to  be  considered  may  be  stated  as  follows.  Given  an 
elastic  potential  W(I),  we  seek  a  suitably  regular  plane  deformation  y 

rv 

on  II,  subject  to  (1.3),  as  well  as  a  nominal  stress  field  o  and  a  pres¬ 
sure  field  p  such  that  (1.12)  holds,  while  £  satisfies  the  equilibrium 
equations  (1.8);  in  addition  y  is  to  be  consistent  with  the  pre-assigned 

/v 

deformation  at  infinity,  whereas  a  should  conform  to  the  boundary  condi- 
tion  for  the  crack  faces.  The  latter,  in  view  of  (1.10),  are  satisfied 
provided 

aa2(xi»Of)  =  0,  oa2(x1>0-)  =  0  (-b <  x^  <  b)  .  (2.2) 

Finally,  we  stipulate  that 


y(x)  •  Fx  +  0(1)  as  x?  +  x?-,oo  ,  det  F  «1 

sw  ^  rw  i  Z  /w 


(2.3) 


where  F  is  a  constant  tensor. 

The  function  y  should  be  twice  continuously  differentiable  on 

/v 

II,  as  required  before,  and  continuous  up  to  £,  while  F  is  to  be  con¬ 
tinuous  up  to  the  interior  of  £.  Further,  we  assume  that  W  is  twice 
continuously  differentiable  on  [2,oo)  and  throughout  its  interval  of 
definition  meets  the  ellipticity  conditions  (1.17).  Thus,  taking  for 
granted  the  existence  of  a  solution  to  the  global  problem  at  hand,  el¬ 
lipticity  is  bound  to  prevail  on  II,  and  £  is  assured  to  have  continuous 
partial  derivatives  of  all  orders  on  II. 

Let  jO  be  the  class  of  all  {y,  a,  p]  that  satisfy  the  field  equa- 

/w  /v 

tions  referred  to  above,  as  well  as  the  boundary  conditions  (2.2).  Clearly, 
a  solution  to  the  complete  problem  formulated  above  belongs  to  £;  con¬ 
versely  any  member  of  &  that  meets  condition  (2.3)  is  a  solution  to  this 
problem.  For  future  purposes  it  is  essential  to  note  that  if  [y,  a,  p] 

/v  /v 

is  in  jS,  the  same  is  true  of  [Qy,Qo,p]  for  every  proper  orthogonal 

rv  <v  /v  /v 

tensor  (^.  This  claim  is  readily  confirmed  because  of  the  objectivity 
of  the  constitutive  law  (1.12)  and  by  virtue  of  the  form  of  the  boundary 
conditions  (2.2). 

At  this  stage  we  recall  certain  properties  of  the  familiar  solutions 
to  the  analogous  crack  problem  within  the  linearized  equilibrium  theory 
of  plane  strain,  which  encompass  both  compressible  and  incompressible  elas¬ 
tic  solids,^-  In  view  of  the  principle  of  superposition  it  is  sufficient 
in  this  connection  to  consider  separately  the  following  two  loading  modes 
at  infinity. 

^See  for  example  Rice  [1]. 


o  o  2  2 

Mode  I:  o^a“o(l),  a22*a  +  o(l)  as  x^  +  x^oo  , 


o  o  o  2  2 

Mode  II:  0n  =  a22  *  °^)  *  ^2*a  +  °(l)  as  x^  +  x^^oo 


(2.4)" 


(2.5) 


Thus  Mode  I  corresponds  to  uniform  uniaxial  tension  perpendicular  to  the 
crack  faces,  while  Mode  II  corresponds  to  simple  shear  parallel  to  the 
crack  faces. 

For  the  Mode  I  loading  the  resulting  displacement  field  is  symmetric 
about  the  x^-axis,  so  that 

u1(x1,x2)  =  {^(x^- x2>,  u2(x1,x2)  *  -u^x^-x^  .  (2.6) 


On  the  other  hand,  the  displacement  field  associated  with  Mode  II  is  anti¬ 


symmetric,  whence 


u1(x1,x2)  = -u1(x1,- x2),  u2(x1,x2)  =  u2(x1,-x2)  . 


(2.7) 


Further,  the  crack  faces  in  Mode  II  undergo  an  in-plane  gliding  deforma¬ 
tion  and  do  not  separate. 

We  now  cite  the  asymptotic  behavior  of  various  field  quantities 
near  the  crack-tips.  To  this  end  let  (r,0)  be  the  local  polar  coordi¬ 
nates  (see  Figure  2)  defined  by 


x^-b*rcos9,  x2*rsin0,  (0Sr<oo,  -it  £0  Sir)  .  (2.8) 


Accordingly,  the  right-hand  crack-tip  corresponds  to  r**0,  while  the 
crack  faces  are  located  at  0  ■  -ir  and  0*it,  respectively. 


Uj,  and  oag  are  the  components  of  the  displacement  and  stress  fields 
in  the  linearized  theory. 
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For  the  Mode  I  loading,  as  r-*0: 


o  11/2  0  0  .  ,  2  0.  ,  ,  1/2.  ^ 

■  —  r  cos-j  (1  -  2v+  sin  y)  +  o(r  )  , 


o  KI  1/2  ,  9  .  2  6]_l  .  1/2. 

u2  ■  —  r  sin  j  [^2(1  -  v)  -  cos  yj+otr  ) 


>  (2.7) 


y 


O  „  -1/2  0  .  0  .  30.^  ,  -1/2.  ^ 

oi:l=  Kj  r  cos-j  (1  -  sin y  sin-^-)  +o(r  )  ,  ^ 


O  „  -1/2  6  30  .  0  ,  ,  -1/2. 

°12“  KI r  cos^cos-sin?+o(r  )  , 


2  wua  2  0x11  2 


V  (2.10) 


o  „  -1/2  0  .  .  0  ,  30.  ,  ,  -1/2. 

°22“  KIr  cos-j  (1+  sin ysin-y)  +  o(r  ) 


y 


Here  y  is  the  shear  modulus  for  infinitesimal  deformations,  v  is 
Poisson's  ratio,  while  is  the  Mode  I  stress-intensity  factor,  given 

by 

KI=a(b/2)1/2  .  (2.11) 


For  the  Mode  II  loading,  as  r-*0: 


O  II  1/2  .  0 

u,  -  -  r  sin-r 

1  V  2 


o  dl  1/2  0 

u2  =  ~T  r  C0S  2 


[2(1-  v)  +  cos2|]  +o(r1/2) 
[-(l-2v)  +  sin2|]  +o(r1/2) 


y  (2.12) 
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o  _  -1/2  .  0  ,  0  0  30  ,  .  ,  -1/2. 

o^j*  Kxx  r  siny  (-2-cosy  cos -y  )+ o(r  ) 

o  „  -1/2  0  ...  .  fl  .  30,,  -  -1/2. 

a12=!  KII  r  cosy  (1  -  siny  sin -y  )  + o(r  ) 

o  _  -1/2  0.0  30  ,  -1/2, 

°22=  Kxx  r  cosy  sin y  cos -y  +  o(r  )  , 


A 


►  (2.13) 

J 


with 

KIi;  =  o(b/2)1/2 


(2.14) 


The  displacements  appropriate  to  the  special  case  of  an  incompressible 
linearly-elastic  body  are  obtained  by  setting  v=l/2  in  the  formulas 
cited  above. 

Although  superposition  no  longer  holds  in  the  nonlinear  theory,  it 
is  still  useful  to  consider  separately  the  kinematic  counterparts  of  the 
Mode  I  and  Mode  II  loadings.  Corresponding  to  the  symmetric  loading 
case  we  now  require  the  deformation  at  infinity  to  be  one  of  uniaxial 
stretching,  characterized  by 


X>  0 


(2.15) 


in  which  X  and  A  ^  are  the  principal  stretch  ratios  at  infinity.  The 
analogue  of  the  antisymmetric  loading  case  is  a  state  of  simple  shear 
governed  by 


(2.16) 
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and  k  is  Che  amount  of  shear  at  Infinity. 

Since  the  governing  field  equations  in  the  finite  theory  do  not 
remain  invariant  if  u^(x^,x2),  u^x^.x^  are  replaced  by  -u^(x^,-x2), 

u2^xl’_x2^'  one  *s  to  won<*er  whether  or  not  the  nonlinear  Mode  II 

problem  admits  a  solution  that  is  antisymmetric  about  the  plane  of  the 
crack.  We  shall  now  show  that,  at  least  for  the  Mooney-Rivlin  material, 
a  solution  with  this  parity  cannot  possibly  exist.  To  this  end  suppose 
now  that  W  is  given  by  (1.19).  Eliminating  £  between  (1.12)  and  (1.8), 
and  invoking  (1.2),  (1.3),  (1.4),  (1.19),  we  find  that 


(2.17) 


Equation  (1.3),  in  view  of  (1.1),  becomes 


1+U1,1+U2,2  +  U1,1U2,2-U1,2U2,1=1  on  11  •  (2'18) 


Now  assume  u  has  the  parity  of  u  in  (2.7).  Then  (2.18)  implies 


1-U1,1~U2,2+U1,1U2,2"U1,2U2,1=1  on  31  '  (2*19) 


Therefore 


Ul,l+U2,2“°  °n  11  ’ 

(2.20) 

i  iuo  9 ~ ui  on  n 

(2.21) 

Mere  and  in  the  sequel  subscripts  preceded  by  a  comma  indicate  partial 
differentiation  with  respect  to  the  appropriate  cartesian  coordinate. 


Next,  P“P  +  P  on  II,  provided 


p'-  ■|[p(x1,x2)-p(x1,-x2)] 


p"-  ■j[p(x1,x2)  +  p(x1,-x2>) 


(2.22) 


By  virtue  of  the  assumed  parity  of  u,  (2.17)  and  (1.1)  give 


1  '  TT 

—  p  =  u  on  II  . 

W  ,a  a,BB 


(2.23) 


From  (2.23),  (2.20)  it  follows  that  p'  is  harmonic  and  u  is  biharmonic 

a 


on  II: 


V2p' » 0  ,  V^u  *  0  on  II  . 
a 

For  the  Mode  II  problem  one  draws  from  (2.16),  (2.3), 

u^»kx2  +  0(l),  u2  =  0(l)  as  p  =  i ^c2  +  x2  -*  co 

We  shall  now  prove  that  (2.20),  (2.21),  (2.24),  (2.25)  imply 
u^«kx2>  u^  ■  c  on  II,  c  ■  constant  , 


(2.24) 


(1.1)  that 


(2.25) 


(2.26) 


so  that  the  entire  deformation  field  is  one  of  homogeneous  simple  shear, 
which  contradicts  the  boundary  condition  (2.2)  for  traction-free  faces 
of  the  crack. 

With  a  view  to  establishing  this  claim,  we  observe  first  that 
u^  -  kx2  and  u2  are  both  biharmonic  on  II,  and  hence  in  the  neighbor¬ 
hood  of  infinity  defined  by  7I"[x|p>b},  2b  being  the  crack  length; 


further  both  of  these  functions  are  bounded  as  p  -* oo .  As  shown  in  the 


Appendix,  u^  then  admits  an  expansion  of  the  form 

00 

u  (x  ,x,)»k6  ,x.+  l  [a^cos  jcp  +  b^sin  jcp+  c^cos(j  -  2)cp 
o  1  2  al  2  j  J  J 

+  dj(a)sin(j  -  2)cp]p_j  on  ,  (2.27) 

where  cp  is  the  polar  angle  defined  by  x^*pcoscp,  X2  =  psincp.  The 

series  in  (2.27)  is  absolutely  and  uniformly  convergent  on  every  closed 

subset  of  71  .  Substitution  from  (2.27)  into  (2.20),  (2.21),  in  view  of 

the  parity  of  u^  ,  confirms  that  these  equations  cannot  hold  unless  all 

(2) 

coefficients  apart  from  a^  in  the  above  series  vanish.  Continuation 

of  the  biharmonic  functions  u  from  71  onto  II  now  yields  the  desired 

a 

conclusion  (2.26). 

In  view  of  the  above  result  it  is  a  plausible  conjecture  that  the 
Mode  II  global  crack  problem  under  consideration  fails  to  admit  an  anti¬ 
symmetric  solution  for  every  legitimate  choice  of  the  plane-strain  elastic 
potential  W. 


1 
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3.  Asymptotic  analysis  of  the  elastostatic  field  near  the  tips  of  the 
crack. 

In  this  section  we  aim  at  the  local  structure  near  one  of  the 
crack-tips  of  the  solution  to  the  global  problem  stated  in  Section  2. 

For  this  purpose  it  is  sufficient  to  consider  the  right-hand  crack-tip, 
situated  at  x^*b,  X2*0.  We  suppose  that  the  material  has  an  elastic 
potential  of  the  power-law  type  (1.18).  Further,  if  (r,0)  are  the  local 
polar  coordinates  introduced  earlier  (see  Figure  2) ,  we  assume  that  the 
deformation  field  of  the  solution  admits  an  asymptotic  representation  of 
the  form 

m  •  i 

^(x^.Xj)  ~r  Ua(0)  as  r-'O  (no  sum  on  a)  .  (3.1) 

Here  m,  and  m2  are  originally  unknown  constants  obeying 

0<m«*  mi^m^.n^]  <  1  ,  (3.2) 

whereas  U  are  as  yet  undetermined  functions  that  fail  to  vanish  iden- 
a 

tically  on  [ — tt , tt ] .  In  addition  we  take  for  granted  the  validity  of  the 

asymptotic  equalities  obtained  by  two  successive  formal  differentiations 

of  (3.1).  Note  that  (3.2)  implies  that  the  deformation  is  continuous  at 

the  crack-tip,  and  that  at  least  one  of  the  deformation  gradient  components 

F  a  becomes  unbounded  there.  Observe  also  that  according  to  (3.1),  (3.2) 
ap 

the  deformation  image  of  the  crack-tip  has  for  convenience  been  placed 


The  asymptotic  equality  symbol  is  used  in  its  standard  connotation. 

Thus,  f(r)~g(r)  abbreviates  f  (r)  -  g(r)  +  o(g(r))  as  r-*0. 
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at  the  origin  of  the  cartesian  coordinates,  which  can  always  be  ar¬ 
ranged  by  means  of  a  suitable  translation.  Finally,  we  suppose  that 
the  pressure  field  associated  with  the  global  solution  satisfies 

p(xlfx2)  ~ri'P(e)  as  r-*0  ,  (3.3) 

Jl  being  another  constant  exponent,  and  assume  that  this  asymptotic 
identity  may  be  formally  differentiated  at  least  once. 

As  pointed  out  in  Section  2,  the  unknown  global  solution  [y,o,p] 
necessarily  belongs  to  the  class  £  introduced  there,  regardless  of  the 
particular  prescription  (2.3)  at  infinity.  It  follows  that  the  determina¬ 
tion  of  y^  and  p  in  accordance  with  (3.1)  and  (3.3)  for  an  arbitrary 
member  of  &  is  bound  to  encompass  the  asymptotic  behavior  of  the  solu¬ 
tion  to  the  complete  crack  problem  at  hand. 

We  now  show  that  for  the  present  purpose  (3.1)  may  —  without  loss 
of  generality  —  be  replaced  by 

y  (x  ,x_)  ~r“u  (0)  as  r-*0  .  (3.4) 

a  i  l  a 

With  this  aim  in  mind  we  first  note  on  the  basis  of  the  discussion  in 
Section  2  that  whenever  [y,a,p}  is  in  &,  the  same  is  true  of 

iv  /v 

[y*,o*,p*}  ,  where 

/V  /Nrf 

y*«  Qy,  cr*mQo,  p*  ■  p  ,  (3.5) 

»v  /v  /v  iv  rv  iv 

and  £  is  an  arbitrary  proper  orthogonal  tensor.  The  components  of 
every  such  tensor  in  the  underlying  coordinate  frame  admit  the  representa¬ 


tion 
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-sln  cp 
cos  cp 


(3.6) 


for  some  cp  in  the  interval  [0,2tt).  From  (3.6),  (3.5),  (3.1)  one 
finds 


m.  m.  .  m  >  .  m  . 

y*  ■  r  U^(9)coscp-r  ^(Ojsin  cp+ o^r  J  +  o^r  J 

“l  m2  /  mU  /  m?\ 

y£  *  r  XU^ (0)sin cp  +  r  ^t^Wcos  cp  +  o(r  j  +  of  r  ) 


>-  (3.7) 


*  y 


whence  there  are  functions  U*  such  that 

a 

y*(x.,x.)*rV(0)  +  o(r”)  .  (3.8) 

a  x  i  a 

Further,  since  is  not  permitted  to  vanish  identically,  it  is  always 

possible  to  choose  cp  so  as  to  assure  that  U*  fails  to  vanish  identi¬ 
cally.  Now,  dropping  the  asterisk  in  (3.8)  we  arrive  at  (3.4). 

Note  that  the  transformation  (3.5)  represents  a  rigid  rotation  of 
the  deformed  body  about  the  origin,  determined  by  the  rotation  tensor  £; 
further,  the  same  rotation  is  applied  to  the  nominal  stress  field  (j  as¬ 
sociated  with  the  original  deformation.  Consequently  the  asymptotic  analy¬ 
sis  to  follow  can  determine  the  local  behavior  of  the  global  solution  at 
most  to  within  an  unknown  rigid  rotation,  which  would  ultimately  depend 
on  the  particular  prescription  at  infinity.  In  this  connection  it  should 
be  emphasized,  however,  that  most  of  the  desired  local  field  quantities 
are  in  fact  invariant  under  such  a  rotation.  For  the  special  case  of  the 
Mode  I  loading  the  local  solution  is  bound  to  be  symmetric  about  the 
x^-axis  and  this  additional  restriction  supplies  essential  information 


l  A 
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relevant  to  the  determination  of  the  rotation  tensor  (£.  As  is  clear 
from  the  discussion  at  the  end  of  Section  2,  analogous  parity  restric¬ 
tions  are  not  available  for  the  Mode  II  loading. 

We  now  seek  to  determine  the  smallest  exponent  m  in  (0,1)  and 
the  functions  Uq  appearing  in  (3.4)  consistent  with  the  governing  field 
equations  and  boundary  conditions.  This  objective  is  most  easily  reached 

by  recourse  to  the  polar  coordinates  (r,0).  Writing  y  (r,0)  in  place 

a 

of  ya(x^(r,0),  X2(r,0))  one  finds  that  the  incompressibility  condition 
(1.3)  becomes 


l/<*l  <*2  <*2  M  . 

r  \  3r  36  3r  3 0  / 


(3.9) 


Further  (1.3),  (1.8),  (1.12)  lead  to 


fc  -  *  (I)Hr+  2S"  (I>  [*rr  f] 

|f  -  2H'  (I)He+  2W"  (X)  [gr„  g09  f§] 


(3.10) 


and  (1.4)  furnishes 


1  “  8rr+  2  800  ’ 


(3.11) 


with 


3y 


H 


a  „2 


8y 


a  r,2. 


7  y„  .  -  as-  7  yn  » 


3y  3y  3y  3y  3y  3y 

q  q  q  q  m  q  q 

8rr  3r  3r  ’  gr0  3r  30  ’  800  ”  30  30 


1 


l 

| 


f 

i 


(3.12) 
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The  boundary  conditions  (2.2)  take  the  form 


o.  3yi  3y2 

2W  (I)  36  +rPTT  =0  ’ 


3y2  3yl 

2W  (I)  — - rp  =0  on  6  *  ±tt 


(3.13) 


With  the  aid  of  (3.12)  and  (3.9),  equations  (3.13)  are  readily  found  to 


be  equivalent  to 


p«  2W'  (D/grr  *  2W'  (I)g00/r2  ,  gr0  =  O  for  0  =  ±ir  . 


(3.14) 


We  thus  need  to  determine  m  and  Uq  in  (3.4)  consistent  with  (3.9), 
(3.10)  and  (3.14). 

Equation  (3.11),  together  with  (3.4),  gives 

I_r2(m-1)j-'2(e)  + •2(0)+m2[u2(0)+u2(e)]]  f  (3.15) 


where  the  dot  denotes  differentiation  with  respect  to  0.  We  shall  hence¬ 
forth  take  for  granted  that  and  U2  do  not  have  a  common  multiple 

zero  on  [  — it , it ]  so  that  the  coefficient  of  r2^m  1)  (3,15)  does  not 

vanish  on 

Next  (3.4),  (3.9)  lead  to 


J-mr2(m“1)[U1(0)U2(0)  -U2(0)U1(0)]  +o(r2(m"1))*  1 


(3.16) 


Dividing  this  identity  by  r2^m  1)^  proceeding  to  the  limit  as  r-*0. 


and  bearing  in  mind  that  m<l,  we  obtain 


L 


-  U2U^  ■  0  on  [-x,ir) 


(3.17) 
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Thus , 


U  =  a  U  on  [-tt.tt]  ,  a  /0  , 

a  a  a 


(3.18) 


in  which  a^  and  are  constants  and  U  an  unknown  function. 

Equations  (3.18),  (3.15),  (3.12),  (3.10),  (1.18)  eventually  lead 


~  2Aa2nr2(m_1)n_1mU(e)Z(e)  , 


f§  ~  2Aa2nr2(m-1^nU(6)Z(0)  , 


(3.19) 


provided 


Z-Gn_2[G(U  +  m2U)+ (n-l)[GU+2m(m-l)GU]]  on  [—tt.it]  , 


„  •2.  2  2  r  .  / 2  .  2 
G*U  +m  U  on  [-tt.ttJ,  a  *  /a^  +  a^ 


On  the  other  hand,  the  boundary  conditions  (3.14)  yield 


(3.20) 


UU  =  0  at  e  =  ±7T  , 


2n  2 (m- 1 ) n  n-1 *  2  . 

p  ~  2Ana  r  G  U  at  0  =  ±ir 


(3.21) 


(3.22) 


We  now  show  on  the  basis  of  (3.22),  (3.21),  (3.19)  that 


Z  »  0  on  [— tt  ,ir] ,  U(-tt)  =  U (tt)  =  0 


For  this  purpose  we  integrate  the  first  of  (3.19)  to  obtain 


_  2n  2(m-i)n  mU(0)Z(0) 

P  ZAa  r  2(m-l)n  ’ 


(3.23) 


(3.24) 
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According  to  (3.21)  either  U(tt)  or  U(ir)  must  vanish.  By  assumption 
U  has  no  multiple  zeros  on  [-rr,n]  so  that  U(ir)  and  U(ir)  cannot 
vanish  together.  Suppose  U(it)  =  0,  U(ir)  #  0.  Then  (3.24)  gives 

p(r,ir)  =»  o(r^m  as  r-»0  ,  (3.25) 

contradicting  (3.22).  Thus  U(ir)  t  0,  U(tt)  =  0  and  similarly  U(-ir)  ?  0, 
U(-tt)  *  0.  On  comparing  (3.24)  and  (3.22)  we  infer  that 

Z(ir)  =  Z(-tt)  =  0  .  (3.26) 

Next,  we  differentiate  (3.24)  with  respect  to  0  and  use  the  second 
of  (3.19)  to  confirm 

UZ  +  jj.  +  — 1~°1— ZU=  0  on  [-Tr,  ttJ  .  (3.27) 

Equations  (3.27)  and  (3.26)  necessitate 

Z  =  0  on  [  —-ft  ,  7T  ]  ,  (3.28) 

which  proves  (3.23).  Thus 

G(U  +  m2U)  +  (n-1)  [GU+ 2m(m-l)GU]  =  0,  U(-ir)  -  U(ir)  «  0  .  (3.29) 

Equations  (3.29)  constitute  a  nonlinear  eigenvalue  problem  that 
arose  also  in  connection  with  asymptotic  analyses  of  other  crack  prob¬ 
lems  in  finite  elasticity  [3],  [4],  [5]  and  was  solved  in  [3].  The 
unique  eigenvalue  in  the  range  (0,1)  is  given  by 


m“l-l/2n,  l/2<n<oo  , 


(3.30) 
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while  the  associated  eigenfunction  U  is  supplied  by 


1/2 

u(e)  -  .10(0/2) [1  -  -12^°^-e(,9n'')2>]  !.(»,.)*< 


k/2 

cos  0]  ,  (3.31) 


where 


2  2 

u(0,n)  -  [1- ic  sin  8]  ,  x=(n-l)/n  .  (3.32) 


The  representation  of  the  local  deformation  field  furnished  by 
(3.4),  (3.18),  (3.31)  yields  only  the  weak  estimates 

J  *  o(r2^m~1^) ,  p»o(r-1)  as  r-0  ,  (3.33) 

and  is  therefore  inadequate,  furthermore,  the  Jacobian  determinant  of 
the  dominant  terms  in  (3.4)  is  found  to  vanish,  which  reflects  the  de¬ 
generate  character  of  the  asymptotic  approximation  established  so  far. 
With  a  view  towards  refining  this  approximation  we  now  replace  (3.4)  by 

y  (r ,9)  =  a  rmU(9)  +  rmV  (0)+o/rm),  m/>m  ,  (3.34) 

with  m*  another  unknown  exponent  and  Vq(9)  unknown  functions  that  are 
not  permitted  to  vanish  identically,  while  m  and  U  are  now  given  by 
(3.30),  (3.31).  We  suppose  that  (3.34)  may  be  formally  differentiated 
at  least  twice.  Again,  by  virtue  of  the  discussion  that  led  to  the  adop¬ 
tion  of  (3.4)  in  place  of  (3.1),  no  generality  is  lost  in  assuming  equal 
exponents  in  the  second  term  of  (3.34). 

From  (3.34),  (3.9)  one  now  draws 
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J  =  r^  ■2(mUy-m/U4')  +  0(riItfm  '2)=.  i  , 


(3.35) 


where 


¥  =  3^2-32^^  on  [— 


(3.36) 


Consequently 


m+m -2s0  , 


(3.37) 


and  further 


mU¥  -  m7  U¥  =  0  on  [— it ,  it  ]  if  m<m7<2-m  , 


(3.38) 


mU¥-m/U¥=‘l  on  [•“*»*]  if  m7»2-m 


(3.39) 


Next,  combining  (3.34),  (3.10),  (3.11),  bearing  in  mj.nd  (1.18),  and  re¬ 
calling  that  U  satisfies  (3.29),  one  confirms  with  the  aid  of  arguments 
similar  to  those  used  in  deriving  (3.28)  that 


Y=0  on  (— 7r ,  tt  ]  if  m<m7<2-m  , 


(3.40) 


Y  * - HU  +  m  U)  on  [—it  ,  tt  ]  if  m7“2-m  ,  (3.41) 

nA 


where  Y  is  an  auxiliary  function  to  be  defined  presently.  In  fact, 
Y*G(x  +  m7  2x)  +  2K(U  +  m2U)  +  (n  -  l)xG 


+  2(n  -  1)  [m7  (m- l)xG  +  m(m7  +  m- 2)UK+UK]  ,  (3.42) 


in  which 

•  •  t 

X*  a^V i  +  a2V ^2»  K-xU  +  o®  XU  on  f  — tt  ,  tt  J 


(3.43) 
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Finally,  the  boundary  conditions  (3.14)  lead  to 

x(±ir)  =  0  if  ,  (3.44) 

¥(±tt)-0  if  m  <  m'  <  2  -  m  ,  (3.45) 

^(ir)  =  l/mU(Tr) ,  'i'(-ir)  *  l/mU(-ir)  if  m,=  2-m  .  (3.46) 

In  view  of  (3.42),  equation  (3.40)  together  with  (3.44)  constitute 
an  eigenvalue  problem  for  x  with  m  as  eigenvalue  parameter.  On  the 
other  hand,  (3.41)  is  an  inhomogeneous  differential  equation  for  x.  to 
be  solved  subject  to  the  boundary  conditions  (3.44).  Consider  next  the 
characterization  of  V  through  the  differential  equations  (3.38),  (3.39) 
together  with  (3.45),  (3.46).  Since  these  boundary  conditions  are  readily 
seen  to  be  implied  by  (3.38),  (3.39),  additional  information  is  needed 
to  determine  the  appropriate  solutions.  Such  information  is  supplied  by 
the  assumption  that  the  solution  V  to  (3.38)  or  (3.39)  possesses  con¬ 
tinuous  derivatives  of  all  orders  on  [  — tt , tt ] .  Equation  (3.38)  now  poses 
an  unconventional  eigenvalue  problem:  one  is  to  determine  the  values  of 
m'  on  the  interval  (m,2-m)  for  which  the  solution  to  (3.38)  is  infi¬ 
nitely  many  times  continuously  differentiable  on  [  —it , tt ] . 

We  turn  at  present  to  the  determination  of  the  relevant  eigenvalues 
of  (3.40),  (3.44)  or  (3.38).  If  at  least  one  of  the  foregoing  problems 
has  eigenvalues  in  the  appropriate  range,  we  choose  for  m'  the  smallest 
such  eigenvalue.  The  following  three  possibilities  now  arise.  In  case 
m'  so  determined  is  an  eigenvalue  for  (3.40),  (3.44)  but  not  for  (3.38), 

X  oust  be  the  eigenfunction  associated  with  m' ,  while  H1  ■  0  on  [  — tt  , tt ] . 


Next,  if  m/  is  an  eigenvalue  for  (3.38)  but  not  for  (3.40),  (3.44), 

¥  has  to  be  the  eigenfunction  associated  with  m/ ,  while  x“0  on 
[-ir,ir].  Further,  if  m/  is  an  eigenvalue  common  to  both  problems, 

¥  and  x  are  the  eigenfunctions  associated  with  this  value  of  m'  . 

Finally,  suppose  neither  of  the  eigenvalue  problems  considered 
above  has  an  eigenvalue  in  (m,2-m).  Then  evidently  m,*2-m  and  the 
inhomogeneous  equations  (3.39)  and  (3.41),  (3.44)  are  to  be  solved  for 
X  and  ¥,  which  in  turn  determine  and  V^. 

At  this  stage  we  investigate  the  solutions  to  (3.38)  and  (3.39). 
Consider  first  (3.38).  In  view  of  (3.31),  we  find 

¥  =  a^U111  on  (0,tt),  ¥*a^(-U)m  on  [-it,0)  ,  (3.47) 

in  which  a^  and  a^  are  constants.  If  ¥  is  to  possess  unlimited 
smoothness  on  [— ir ,tt]  »  m‘  /m  must  be  a  positive  integer.  Thus  the 
smallest  eigenvalue  m' >m  is  given  by 

m/«2m*2-l/n  ,  (3.48) 

and  the  associated  eigenfunction  by 

¥-bQU2  ,  (3.49) 

where  bg ^  0  is  a  constant.  In  addition  n/  defined  by  (3.48)  conforms 
to  (3.37),  provided 

n <  3/2  .  (3.50) 


In  the  event  that  n2  3/2,  (3.38)  does  not  possess  an  eigenvalue  in 
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(m,2-m).  With  a  view  to  the  determination  of  ¥  when  ni3/2  we  con¬ 
sider  (3.39)  merely  on  the  interval  [0,it],  which  will  presently  be  seen 
to  be  sufficient.  Evidently 

on  [-ir.ir]  ,  (3.51) 

if  one  sets 

V1  (9)  “'jt'Ke)  -  t(-0)  ] ,  H'//(0)*^['i,(0)+  4'(-6)  ] ,  -it  £  0  £  it  .(3.52) 

Moreover,  v'  and  v"  are  found  to  satisfy 

mU^ - m' U'?'  -  0  on  [O.tt]  ,  (3.53) 

mU1?" -  m' UY"  -  1  on  [0,it]  .  (3. 54) 

If  ni3/2,  the  only  solution  to  (3.53)  having  an  infinitely  smooth  odd 
extension  to  [-n,ir]  is 

y'=0  on  .  (3.55) 


Further,  (3.54)  gives 


V"  (8) 


i  r _ deg _ 

m)„  [U(<p)]2/m 


m'  =  2  -  m,  (0  £  0  £  it)  , 


(3.56) 


with  b^  an  arbitrary  constant.  We  now  change  the  variable  of  integra¬ 
tion  in  (3.56)  from  <p  to  £  by  means  of  the  one-to-one  mapping 

2  1/2 
1  [l  +  tcsin  (p  -  a>(tp.n)cos  cp] 

n/2  utCcp* n)  +  tccos  ^ 


cos  C 


(0scsf)  .  0.57) 


u)  and  k  being  given  by  (3.32).  This  change  of  variable  carries 
(3.56)  into 


Y"  (0)  -  U1 


V  /tn 


bl- 


.  1-1/m 


2nT 


f^0  1 + 2(n-l)cos2g 


[cos  Cl 


2/m 


dS 


(3.58) 


where 


cos  5. 


2  1/2 
1  [1+Ksin  9  -  a)(6.n)cos  9] 


n/2 


j(0.n) +  fCCOS  0 


(0  S0  Sir) 


(3.59) 


The  integral  in  (3.58)  may  be  evaluated  in  terms  of  hypergeometric  functions, 
provided  n>3/2.  Following  an  analysis,  the  details  of  which  will  be  omit¬ 
ted  here,  one  is  led  to  the  representation 


b2Um/  /m  +  H 


on  [0,ir] , 


m 


2  -  m 


(3.60) 


Here 


.  4:n3''1/m  r  ( 1/ 2  -  1/m) 

62"Dl"  2m2  r(l-l/m)  * 

r  denoting  the  Gamma  function,  while 

5/2 

H(9)  =  -  n  ~2  [u(0,n)  +  kcos  0] 


(3.61) 


X  F(l/2  -  l/m,l/2;3/2  -  l/m;cos2Cg)  -  t«sin  j  .  (3.62) 


i  a 
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in  which  F  stands  for  the  hypergeometric  function. ^  Now  if  v"  is 

to  possess  continuous  derivatives  of  all  orders  on  fO.ir]  then 
2 

must  vanish.  From  (3.51),  (3.55),  (3.60)  one  infers 

'f'(e)  *  H(0)  (-Trsesr)  .  (3.63) 

One  draws  from  (3.62),  (3.63)  and  familiar  properties  of  the  hypergeo¬ 
metric  function  that 


¥(0)  =-2n1/2(2mrm/m'  , 


*(±ir) 


/mr  nm  T (3/2  -  1/m) 
nm/  ^(1 ~  1/m) 


►  (3.64) 


the  first  of  which  is  also  a  direct  consequence  of  the  differential 
equation  (3.39).  The  second  of  (3.64)  will  be  used  later  in  discussing 
the  shape  of  the  deformed  crack. 

In  case  n™  3/2,  (3.39)  fails  to  admit  a  solution  of  the  requisite 
smoothness.  We  postpone  consideration  of  this  special  case,  which  requires 
separate  attention,  until  Section  4  and  turn  to  the  determination  of  x* 
Consider  first  (3.40)  and  set 


X“X'  +  x"  on  [-n,n] 


(3.65) 


where 

X  (9)  *-|[x(0)  -  X  C— 8)  1 ,  x"  (9)  ■  jfxO)  +  x(-9)  J  (-use  in)  .  (3.66) 


1 


See  for  example  [15],  p.556. 


^Note  that  H(0) 
of  all  orders  on 
0  -  0. 


is  an  even  function  of  0 
[-n,n]  and  recall  that 


with  continuous  derivatives 
U(0)  has  a  simple  zero  at 
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From  (3.65),  (3.66),  (3.40),  (3.42)  we  have,  bearing  in  mind  the  parity 
of  U, 

G(x'  +  m'  V)  +  2K'  (U  +  m2U)  +  (n-l)x'G  +  2(n-l)  [m'  (m  -  l)x'G 

+  m(m'  +  m-  2)UK/  +  UK'  ]  -  0  on  [0,ir]  ,  (3.67) 

G(x//+m/2x//)  +  2K"  (U  +  m2U)  +  (n  -  1)  x"G+  2(n  -  1)  [in' (m  -  1)  x"G 

+  m(m'+m  -  2)  UK"  +  UK"  ]  -  0  on  [0,ir]  ,  (3.68) 

where 

K'- x'u  +  mm'x'u,  k"  =  x" U  +  mm'  x" U  .  (3.69) 

Next,  (3.66),  (3.44)  furnish 

x'(*)-0  ,  x'(0)«0  ,  (3.70) 

x"  (iO-O  ,  x"  (0)*  0  .  (3.71) 

The  differential  equation  (3.67)  and  the  boundary  conditions  (3.70) 
are  also  encountered  in  [4].  As  indicated  there,  the  transformations 


1  [1  +  Ksin29  -  u)(8  ,n)cos  8  ]  ^2 

n/J  u)(9,n)  +  tccos  0 


(OS0  £  it) 


(3.72) 


/ 

W^(£)  *  [w(8,n)  +  iccos  9]  ”  x#  <«)  (0  s  0  Str)  ,  (3.73) 

reduce  (3.67),  (3.70)  to 

WX(C)  +  X2W1(C)  -  0  (0*C*"/2),  W1(0)  -  W^ir/2)  -  0 


(3.74) 
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where 

A^  *  4nmy  (nm*  -  n  +  l)/(2n  -  1) 

The  transformed  eigenvalue  problem  (3.74)  has  the  solutions 
WX(C)  -  t^cos  A  £,  A-  j  (j  -  1,3,5,..,) 

in  which  is  a  constant.  Under  the  transformation 

/ 

W2(£)  =  [u>(0,n)  +  xcos  8]-m  x”  (9)  (0^0  sir)  , 

together  with  (3.72),  equations  (3.68),  (3.71)  become 

W2(C)  +  A2^^)  -  0  (0S£  Sir/2),  W2(0)  =W2(ir/2)  «  0  , 

and  the  solution  to  (3.78)  is  given  by 

W2(0  -  c2cos A £ ,  A  =  j  (j  *  2,4,6, .. .)  . 

Thus  the  eigenvalues  of  the  original  problem  (3.40),  (3.44) 
the  solutions  for  m  of 

4nm/  (run'  -  n  +  1)  /  (2n  -  1)  *  (j*  1,2,3,...) 

For  j  an  odd  positive  integer,  x"  ”  0  on  [-ir.ir]  and  x*(9) 
nished  by  (3.76),  (3.73),  (3.72),  while  if  j  is  even,  x#(9)  =  0 
[-ir,ir]  and  x"  (0)  is  furnished  by  (3.79),  (3.77),  (3.72). 

On  discarding  negative  roots  m  of  (3.80)  one  arrives  at 


(3.75) 

(3.76) 

(3.77) 

(3.78) 

(3.79) 

are 

(3.80) 

is  fur- 
on 
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Fur  ther, 

m'«=l-l/2n*m  if  j  -  1  ,  (3.82) 

so  that  j  2-  2  since  m^m  .  The  minimum  admissible  eigenvalue  xa 
occurs  at  j  »  2  and  is  given  by 

m/ * [ic  +  VK+8m/n  ]  ,  K=(n-l)/n  ;  (3.83) 


the  corresponding  eigenfunction  is  readily  shown  to  be 


X(0)  ”  c2[a)(0,n)  +  kcos  0]m 


1+Ksin  6  -  to (Q  ,n)cos  0 

2  2  ~  1 

[oi(0,n)  +  kcos  0]  n 


(-it  ^  0  ^tt)  .  (3.84) 


The  eigenvalue  (3.83)  conforms  to  (3.37)  as  long  as 


n<  7/2 


(3.85) 


Finally  suppose  that  n2  7/2,  in  which  case  (3.40),  (3.44)  do  not 
possess  an  eigenvalue  on  (m,2-m).  Equations  (3.41),  (3.44),  (3.65), 
(3.66)  in  conjunction  with  the  transformations  (3.72),  (3.73),  (3.77) 
lead  to 


W1(5)  +  ^2W1(C) 


cos  £  ,  W^(0)  =Wj(ir/2)  *  0 

m  A 


w2(5)  +  \2v2  (C)  -  0,  W2(0)  -  W2(tt/2)  -  0  , 


>  (3.86) 


J 


where 


X2-3(2n+l)/(2n-l),  K«(n-l)/n  . 


(3.87) 
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Thus, 


wx(0  = 


-Bicni^n 

4Am^(l+n) 


cos  £ 


if  n  a  7/2 


(3.88) 


W2(£) *  c2cos 2  £  ^  n*7/2  (c2  *  constant)  , 

W2(C)«0  if  n>  7/2  . 

Equations  (3.88),  (3.89),  together  with  (3.77),  (3.73),  (3.72), 
furnish 


►  (3.89) 

J 

in  turn 


Bic/In 


X(9)  * - "9*  —  [to ( 0  ,n)  +  kcos  0]m  1  /l  +  <sin^0  -  a>(0  ,n)cos  0 

4Am  (1+n) 


+  CjfuO.n)  +  tccos  6] 


m 


1+Ksin  9-a>(9,n)cos9  . 

2  2-i 

[u>(9 ,n)+KCOS 0 ]  n 


(-it  S0  Sir), 
m  *  2  -m  , 


(3.90) 


for  ni7/2,  with  c2«  0  when  n>7/2. 

At  this  point  we  summarize  the  various  possibilities  that  arise 
in  the  determination  of  the  second  term  in  (3.34).  First,  if  nz7/2, 
there  are  no  eigenvalues  of  (3.38)  or  (3.40),  (3.44)  on  the  interval 
(m,2-m).  Thus  m/*2-m  in  this  instance  and  x*  'f  are  given  by 
(3.90),  (3.63).  On  solving  (3.36)  and  the  first  of  (3.43)  for  Vq,  one 
is  led  to 


V 


llx  “  a2*^a 


V2-  (a^  +  a2x)/a 


on  [  —IT  ,  n  ] 


(3.91) 
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Next,  suppose  3/2<n<7/2.  In  this  case  only  the  x~problem  admits 
an  eigenvalue  on  (m,2-m),  which  is  given  by  (3.83),  whereas  x  is  sup¬ 
plied  by  (3.84).  is  now  found  by  setting  ¥  ■  0  in  (3.91),  which 

gives 

Vi«  aiX/s2,  V2=a2x/a2  on  [-rr.ir]  .  (3.92) 


Finally,  for  the  range  l/2<n<3/2,  in  which  case  both  the  x~prot>lem 
and  the  ^-problem  admit  eigenvalues  on  the  appropriate  range,  we  seek 
the  smallest  such  eigenvalue.  This  criterion,  along  with  (3.48),  (3.83), 
leads  to 


m'  =  +  8m/n j  if  l<n<3/2 

m7  -  1  if  n  *  1  , 


>  (3.93) 


m'»2-l/n  if  l/2<n<l 


J 


Further, 


¥  *  0  on  [-ir,ir]  for  l<n<3/2,  x*0  on  [—it  ,  tt  ]  for  l/2<n<l  .  (3.94) 

For  l/2<n<l,  Va  is  supplied  by 

V1»-a2'y/a2,  V2=a^4'/a2  on  [-tt.tt]  ,  (3.95) 


where  ¥  is  given  by  (3.49).  When  n=l,  both  problems  have  the  common 
eigenvalue  m' *1.  In  the  following  section  we  shall  explore  in  greater 
detail  the  Mooney-Rivlin  material,  which  is  a  power-law  material  with 
n  ■  1. 


It  is  convenient,  for  reasons  that  will  be  made  clear  presently, 
to  apply  an  additional  rigid  rotation  to  the  deformation  field  charac¬ 
terized  by  (3.34).  In  particular,  consider  the  rotation  tensor  whose 


component  matrix  in  the  underlying  frame  is  given  by 


a 


+  a 


2 

2  * 


(3.96) 


Calculating  y*  in  the  first  of  (3.5)  from  (3.34),  (3.96),  bearing  in 

<v 

mind  (3.36),  (3.43)  and  dropping  the  asterisk  from  the  resulting  estimates 
for  the  components  y*  ,  one  arrives  at 


yi(r,0)--irmne)+o(rm) 


y2(r,0)  =  armU(0) +-^rm  x(6)+o(rm) 


V  (3.97) 


The  significance  of  (3.97)  is  as  follows.  Suppose,  in  particular,  that 
n>7/2  and  consider  the  deformation- image  of  the  extended  crack  axis 
0-0.  From  (3.97),  (3.90),  (3.64),  (3.31)  follows 


o  l/2._  x-m  /  /  _/ 

y1  (r ,0)  *  — — 7^-  m-  rm+o(rm),  y„(r,0)=o(r  )  , 

J*  in  a  ^ 


(3.98) 


m  a 


and  thus  the  ray  0-0  after  deformation  is  tangent  to  the  x^-axis  at 
the  origin. 

In  case  l<n<7/2  (n^3/2),  ¥  -  0  on  [-ir,ir]  and  the  first  of 

/ 

(3.97)  yields  merely  the  weak  estimate  y^-o(r  ).  There  is,  however, 
no  difficulty  in  establishing  a  non-degenerate  estimate  for  y^.  Indeed, 
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let 


y  (r,0)  =  -  —  rm  X(0)  +  o(rm  ),  m">m  , 
la 


(3.99) 


and  suppose  that  is  given  by  the  second  of  (3.97).  Entering  (3.9) 

with  this  Ansatz,  one  has 


mXU  -  m,/  XU  =  0  if  m<m//<2-m  , 


mXU  -  m"  XU  =  1  if  m"  =  2  -  m  . 


(3.100) 


These  equations  are  simply  (3.38),  (3.39)  with  m'  ,  T  replaced  by  m"  , 
X,  and  hence 


m,/=2--^,  X=bQU2  on  [-■ tt.tt]  for  l<n<3/2  , 


to"  =  1+-^-,  X=H  on  [-it, it]  for  n>3/2 
zn 


Thus  for  l/2<n<3/2. 


(3.101) 


y1(r,0)  = b0r2m[U(6)]2+  o(r2m)  , 


y.(r,0) -armU(0)+-rm/x(0)+o(rm/)  , 

i.  a 


y  (3.102) 


while  for  3/2<n<7/2, 


y. (r,0)  *  -  —  r2-mH(0)  +  o(r^-m  )  , 

X  d 


2-m  < 


y„(r,0)  “  armU(0)  +  —  rm  x(0)+o(rm)  , 

t  a 


y  (3.103) 


m' ,  x  and  H  being  given  by  (3.83),  (3.84),  (3.62).  In  case  n2  7/2, 
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y^(r,0)  =  -  ~r^  -  mH(0)  +  o(r2  -m)  , 

y2(r,0)  =  ar^KO)  +^r2  ~mx(0)  +  o(r2  m)  , 


> (3.104) 


J 


with  x  now  furnished  by  (3.90).  U  and  m  are  supplied  in  all  cases 
by  (3.30),  (3.31).  As  mentioned  earlier,  the  special  case  n=3/2  will 
be  dealt  with  in  the  next  section. 

We  now  turn  to  the  determination  of  the  local  structure  of  the  pres¬ 
sure  field,  which  has  been  assumed  to  admit  the  representation  (3.3).  From 
(3.9),  (3.10),  (3.12)  follows 


2W'(I)72y1+  2”"  W  (  3^ 


+  1 


31  _____ 

3r  T  2  30  30  )  r \3r  30  30  3r 


(3.105) 


Suppose  first  l/2<n<3/2,  in  which  case  y  is  given  by  (3.102).  Sub- 

a 

stituting  from  (3.102)  and  (3.3)  into  (3.105),  one  is  led  to 


/  . 2n-2„n  -1,  .  -1.  nri-£.-2/n  •  •  .  ,  ,  m+t-2.  , 

4nAbga  G  r  +o(r  )  =>  ar  (iPU -mPU)  +  o(r  )  ,  (3.106) 


where  G  is  given  by  (3.20).  By  virtue  of  the  boundary  conditions 
(3.14),  along  with  (3.102), 

p(r,±ir)  ~2nA[maU(*)]2(n~2)r(2-n)/n  ,  (3.107) 

\ 

whereas  from  (3.3), 


Since  bg  0  one  has 


p(r,±ir)  ~r*’P(±ir) 


(3.108) 
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££l-m=l/2n  . 

From  (3.107),  (3.108),  (3.109)  one  now  draws 

P(±w)-0  . 

Suppose  £<l/2n.  Then  (3.106)  gives 

£PU  -  mPU  =  0  on  [-u,ir]  , 


which  together  with  (3.110)  yields 

P  ■  0  on  [-ir,ir] 


Therefore 


£  =  l/2n  , 


(3.109) 


(3.110) 


(3.111) 


(3.112) 


(3.113) 


£PU  -  mPU  ■  4nAbQa2n"3Gn  on  [-ir,ir] ,  P(±ir)  =  0  .  (3.114) 

The  solution  to  (3.114)  is  found  to  be  given  by 

P(0)*4bgAa2n  3n3^2  nm2n  2(u)(0,n)  +  kcos  0]^2n  /l  -  cos2^q  ,  (3.115) 

with  cos  ?0  and  ui(0,n)  given  by  (3.59),  (3.32). 

When  n>3/2,  substitution  from  (3.103)  or  (3.104)  into  (3.105) 
eventually  yields 


£ 


1 


(3.116) 


mPU  -  £PU ■ 


„  ,  2n-4_n-2.  r  , 

2nAa  G  S  on  [-iv,tt] 


(3.117) 
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where 


S  -  G[H  +  (2  -  m) ^3  +  (n-  1)  [GH  +  2(2  -  m)  (m  -  1)GH] 
The  boundary  conditions  (3.107),  (3.108)  give 

P(±ir)  “  2nA(amU(ir)  ]  ^ 


(3.118) 


(3.119) 


The  boundary-value  problem  consisting  of  (3.117),  (3.119)  fails  to  admit 
a  solution  with  continuous  derivatives  of  all  orders  at  0=»O,  and  con¬ 
sequently  the  Ansatz  (3.3)  cannot  possibly  be  consistent  with  the  field 
equations  when  n>3/2. 

In  summary,  the  behavior  of  the  pressure  near  the  crack-tip  is 
governed  by 


p(r  ,0)  =  r^/,^nP(0)  +  o(r^^n)  for  l/2<n<3/2  , 

p(r,0)«o(r  ^+1/n)  f0r  n>3/2  , 


>  (3.120) 


J 


with  P  furnished  by  (3.15);  the  weak  estimate  for  p  appropriate  to 
n>3/2  follows  from  (3.10)  together  with  (3.103)  or  (3.104).  A  non¬ 
degenerate  pressure  estimate  valid  for  n>3/2  would  presumably  neces¬ 
sitate  a  higher-order  asymptotic  analysis  of  the  deformation  field. 

We  note  that  if  n>3/2  the  dominant  terms  in  (3.103),  (3.104) 
have  a  Jacobian  determinant  J~1  and  hence  constitute  a  locally  one-to- 
one  mapping.  In  contrast,  J  computed  from  the  dominant  terms  in  (3.102), 
which  apply  to  l/2<n<3/2,  is  merely  o(r^  ^^n)  .  The  discussion  in 
the  next  section  of  the  Mooney-Rivlin  material  (n*l)  and  of  the  special 
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case  n*3/2  will,  however,  lead  to  an  improved  approximation  to  the 
deformation  at  least  when  lsn<3/2. 


4.  The  transition  case  n*3/2.  Higher-order  analysis  for  the  Mooney- 
Rivlin  material. 

We  first  investigate  the  modifications  necessary  in  the  asymptotic 
structure  of  the  deformation  and  pressure  fields  when  n=3/2.  The  fol¬ 
lowing  analysis  parallels  that  in  [4]  for  the  derivation  of  the  logarithmic 
term  in  the  deformation  field  when  the  hardening  parameter  there  takes  the 
value  7/6. 

Recall  from  the  first  of  (3.103)  that 

y.Cr.e)  --ir2_mH(e,n)  +  o(r2"m),  (n>3/2)  ,  (4.1)1 

X  3 

and  note  on  the  basis  of  (3.62)  that 


H(0,n) 


cQ[U(0,3/2)]' 
n-  3/2 


+  0(1) 


n-y  , 


(4.2) 


Thus  y^(r,9)-*oo  as  n-»3/2  from  above  for  all  0^0  and  every  suf¬ 
ficiently  small  r>0.  Consequently  the  preceding  estimate  for  y^  is 
inadequate  near  n  *  3/ 2 . 

We  shall  now  show  that  in  a  deleted  neighborhood  of  n*3/2, 

y1(r,9)  -  -  jr1+1'/2nH(0,n)  +  bQU2(e,n)r2  1/,n  +  o(r*)  ,  (4.3) 


with  m#  ■  max(l +  l/2n,  2  -  1/n)  and  b^  i  0  a  constant.  Suppose  first  that 
n>3/2  and  replace  (4.1)  by 


Hfe  now  write  U(0,n),  H(9,n)  in  place  of  U(9),  H(0)  in  order  to  em¬ 
phasize  the  dependence  of  these  functions  upon  the  hardening  parameter. 
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y1(r,0)--ir2'1IIH(e,n)  +  r“  T(6,n)+o(rm  ),  m">2-m  .  (4.4) 


Substitution  from  (4.4)  and  the  second  of  (3.103)  into  (3.9)  leads  to 

mUT-m"TU*0  on  [-ir,n]  ,  (4.5) 

which  is  merely  (3.38)  with  m7  ,  41  replaced  by  m"  ,  T.  Thus 

m"  =  2m  =  2  -  1/n,  T  =  bgU2  .  (4.6) 

In  case  l<n<3/2,  a  similar  procedure  yields 

y1(r,0) -bQr2mU2(0,n) --^r2'mH(6,n) +o(r2_m)  ,  (4.7) 


which  confirms  (4.3).  It  should  be  noted  that  the  first  term  of  (4.3) 
dominates  as  r-*0  for  n>3/2,  while  for  n<3/2  the  second  term  dom¬ 
inates;  in  either  case  the  respective  dominant  term  agrees  with  our  pre¬ 
vious  results.  The  two  terms  trade  the  dominant  role  as  n  passes 
through  3/2,  the  respective  powers  of  r  having  the  common  exponent 
4/3  at  n  *  3/2. 

From  (4.2),  (4.3)  now  follows 


y1(r,e) 


1+1 /2n 

c  r 
0 

n  -  3/2 


+  b 


2  -  1/n 

0 


+  rl+l/2n 


n-^e) 


+  bQ(n- 3/2)r2_1/nn2(0) +  o(r”*)  as  n-3/2  .  (4.8) 

where  n^O)  and  are  functions  that  stay  bounded  as  n-*3/2. 

Thus  if  y^(r,0)  is  to  remain  bounded  as  n-*3/2  one  must  have 
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c 

b0*a  7^172  +0(1)  as  •  (4.9) 

Proceeding  to  the  limit  in  (4.8)  as  n-*3/2  one  arrives  at 

y.(r,e)--^  [U(6,3/2)]2r4/3logr  +  r4/3n(0)+o(r4/3)  ,  (4.10) 

2an 

in  which  n(0)  is  bounded.  This  result  motivates  the  following  modifi¬ 
cation  of  (4.3)  in  case  n«  3/2.  We  suppose 

y1(r ,6)  -  r^3log  r  v^S)  +  r^3v2(8)  +  o(r4^3)  ,  (4.11) 

and,  recalling  (3.102),  (3.103),  assume 

y2(r,0)  -ar2/3U(0,3/2)+O(r2/3+6)  ,  (4.12)1 

where  6>0.  Substitution  from  (4.11),  (4.12)  into  (3.9)  yields 

Uv1~2Uv1*0,  Uv2  -  2Uv2  »  2^j.O  ~ on  [-*♦"]  •  (4.13) 

The  first  of  (4.13)  gives 

v^”c^U^  on  [— tt  ,  tt  ]  ,  (4.14) 

so  that  the  second  of  (4.13)  may  now  be  written  as 

Uv2  -  2Uv2  ■ -|^UU^  -  ^  on  [-ir,ir]  .  (4.15) 

Hlerely  requiring  the  error  to  be  o(r^^)  is  inadequate  for  the  present 
purpose. 
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If  (4.15)  is  to  possess  a  solution  of  unrestricted  smoothness  on 

2 

[-w,ir],  then  c^-3cg/2an  in  agreement  with  (4.10). 

Next,  observe  that  In  view  of  (4.9),  (3.115),  P  becomes  unbounded 
as  n-*3/2  for  — it<0<it.  The  foregoing  analysis  suggests  that  we  re¬ 
place  (3.3)  by 


p(r,0)  ~r*/^log  r  P(0)  for  n-3/2  .  (4.16) 

Substitution  from  (4.11),  (4.12),  (4.16)  into  (3.105)  gives 

i  .  2  •  3/2 

-|PU--|PU»6Ac0G  '  on  [-ir,ir]  .  (4.17) 

Moreover,  from  the  boundary  conditions  (3.14), 

P(+ir)-0  .  .  (4.18) 


Consequently , 

P(0)  ■  4cqA[u(0,3/2) +-jcos  J\  -  cos^Cq  for  n*3/2  ,  (4.19) 

2 

in  which  cos  Sq  is  given  by  (3.59)  with  n=>3/2. 

Equations  (4.16),  (4.19)  supply  a  non-degenerate  estimate  for  the 
pressure  in  case  n-3/2.  In  addition,  (4.3)  together  with  the  second  of 
(3.102)  yield  a  locally  one-to-one  approximation  of  the  deformation  if 
1  < n<  3/2. 

We  turn  now  to  the  Mooney-Rivlin  material  (n-1),  which  has  the 
complete  elastic  potential  (1.19),  and  seek  improved  estimates  for  the 
corresponding  deformation  and  pressure  fields. 

The  equilibrium  equations  (3.10)  in  this  instance  reduce  to 


with  H  and  H  supplied  by  (3.12).  The  boundary  conditions  (3.14), 
r  0 

in  turn,  at  present  become 


grQ  “  0  at  9  -  ±ir 


(4.21) 


Setting  n=l  in  (3.30),  (3.31)  one  finds  that 

m-  1/2,  U(0)  *  sin(0/2)  .  (4.22) 

We  recall  also  that  the  eigenvalue  problems  determining  x»  f  —  and  hence 
in  (3.34)  —  admit  the  common  eigenvalue  m'-l  .  Upon  substituting 
n»l  in  the  corresponding  solutions  (3.49),  (3.84)  one  arrives  at 

X=*b^cos0,  ’P*b2Sin^-j  ,  (4.23) 

where  b^  and  b2  are  constants.  Finally,  (3.113),  (3.115)  supply  a 
first  estimate  for  the  pressure: 


*  -  1/2,  P (6) 


2pb2 

a 


0 

cos  y 


(4.24) 


Equations  (3.3),  (3.34),  (3.91)  together  with  (4.22),  (4.23),  (4.24) 
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1/2  6  1  2  0 
y1(r,0)*a^r  '  sin-|  +  — 2  ta^b^cos  0  -  a2b2sin  ^-]r  +  o(r) 

a 


y *  a2r^^sin-|  +  -^-  [a^^sin2-!^  ajb^cos  0]r  +  o(r)  ,  > (4.25) 


p<rj6,.--^ti«co99  +  o(ri/2)  _ 

a 


With  a  view  to  refining  these  estimates  we  first  replace  (3.34)  by 


y  (r ,0)  *  a  r1/2U(6)  +  rV  (0)  +  r”*R  (0)  +  rm*S  (0)  +  o(rm*) 
a  at  a  a  a 


(4.26) 


where  ml>m.  >1  and  R  ,  S  are  as  yet  undetermined,  whereas  the  func- 
"  ”  01  cx 

tions  U,  are  already  known.  Combining  (4.26)  with  (3.9),  (4.20)  and 
invoking  the  boundary  conditions  (4.21)  one  finds  after  considerable  com¬ 


putations  that 


*  3/2,  m^  *  2  , 


(4.27) 


R-L*-^-[a1x1  -  a2'1'1l »  r2  °‘_T[a2X1  +  on  - 28) 


*'^-fa1X2  ~  a2V2^  ’  S2  “_Va2x2  +  al'f,2^  on  ♦  (^-29) 


with  xx.  V2’  x2  suPPlied  by 


1 


_  .  3  0  2blb2  .  0  .  .  2  0  0  4  3  0 

^(0)  "  c2sin  J+ 2 — sin-^-4sin  ycos-^-^cos  j  , 

“a 


.  30  2  .0 

X.(0)  *  c^sin-^ - 2slnT  ’ 

i  i.  2a 


(4.30) 
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f2(e)  -  dlSin4|  +  ^-  (3b1c2 - |+— 6b2c1')sin2| 

a  \  a  a  V 


bib2  8b!  0  3  0  2bl  3  0  6 

+  — £ - ^ cos  y  sin  y - 2* cos  jsiny  , 

a  a  a 


x2^0)  =  — f^sin  0  -  sin  20]  +  d2cos 


20 


KW 

b2 

7  blb2  9 

2  '  +  4 

cos  0  “  2 

4  2  +16  2 

a 

a 

a 

M4.31) 


Here  ,  c2>  d^,  d2  are  constants. 

In  order  to  calculate  an  additional  term  for  the  pressure  field, 
we  apply  the  rotation  with  the  component  matrix  (3.96)  to  the  deformation 
field  (4.26).^  This  yields 


y^(r ,0)  ■  -  ^ry(e)  r3/2'i>1(d)  --^r2'l'2(0)  +  o(r2)  , 


y2(r,0)  -  ar^lK©)  +  -^rx(6)  +-^r3^2x1(0)  r2x2(9)  +o(r2) 

Next,  for  the  Mooney-Rivlin  material,  (3.105)  reduces  to 


>(4.32) 


uV2  .1 

M  yl  r  \ 3r  36  30  3r  j 


(4.33) 


We  now  suppose  that  the  pressure  field  conforms  to 


Kecall  that  the  pressure  field  is  not  affected  by  such  a  rotation. 
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p(r,9)  *rilP(e)  +  rJi  +  ),  l‘ >  i  ,  (4.34) 


with  i,  P  supplied  by  (4.24)  and  l'  ,  P^  unknown.  Equations  (4.32), 
(4.33),  (4.34)  and  the  boundary  conditions  lead  to 


i'-l,  P1(0)=^ 


(4.35) 


Finally,  suppose  that  the  Mode  I  loading  (2.15)  prevails  at  infinity, 
in  which  case  the  local  deformation  field  (4.26)  is  to  be  symmetric  about 
the  x^-axis.  In  this  case, 


a.^0,  a-*  a,  b.  *  0,  c,  «  0,  d„  =  0  , 


(4.36) 


and  (4.26)  reduces  to 


,  1  .  .  20xl  3/2 

y1(r’0)*-7rb2sin  2+^-r 


.  ,  2  0  0  4  3  0 

4sin  y  cosy +  -j  cos  y 


1  2 
•—  r 
a 


4  0  2b2  /.  b2  \  .  2  0 

Vin  2+-j(3cl-^2jsin  2 


+  o(r  )  , 


,  ^  1/2  9.1  3/2 

y2(r,0)*ar  siny  +  yr 


.  30  b2  .  0 

ClSinT-72S±n2 

2a 


+  -y  r3[2sin  0  -  sin  20]  +  o(r^) 
aJ 


L.  (4.37) 


Further,  the  associated  pressure  field  conforms  to 
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2\A>2  - 

p(r,0)  * - 2~ cos ~2  +  — j r  (3  -  c°s  6)  +  o(r)  .  (4.38) 

a  a 

The  above  results  for  the  symmetric  loading  mode  are  identical  with  those 
arrived  at  in  an  unpublished  study  by  Knowles  and  Sternberg,  which  is  sum¬ 
marized  in  [9].  In  contrast,  it  is  not  possible  to  specialize  the  con¬ 
stants  in  (4.26)  so  that  the  displacement  field  u  =y  -x  has  the  parity 

a  a  a 

of  u  in  (2.7). 
a 

This  fact  is  consistent  with  the  conclusion  reached  at  the  end  of 
Section  2,  according  to  which  the  Mode  II  problem  for  the  Mooney-Rivlin 
material  fails  to  admit  an  antisymmetric  global  solution. 
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5.  Discussion  of  the  deformation  and  stresses  near  the  crack-tip. 

Our  initial  objective  here  is  to  explore  the  asymptotic  results 
for  the  deformation  field  established  in  the  preceeding  two  sections.  Of 
particular  concern  is  the  deformation  image  of  the  crack  faces  0  =±ir  . 
Suppose  first  that  n>3/2  and  consider  (3.103).  Then,  to  leading  order, 


y^r.e)  *  -^1+1/2nH(e)  +  o(r1+1/2n)  , 


,(r,6)  =arX  1/2nU(0)  +  o(r1-1/2n)  . 


A 


>  (5.1)' 


In  view  of  (3.31),  (3.64),  one  has  in  particular 


..  /-  .  x  _  1  _l+l/2n  ^  r( 3/2-1/m) 

yi(r,±V)=--r  2"'"  l'/2n  'r(2"-l/m) 

2m  m  n 


,  ,  1+1 /2n.  .  1 

+  o(r  ),  m-  1- ^  , 


,  s  l-l/2n  -tc/2  ,  l-l/2n. 

y2(r,ir)=»ar  n  +o(r  )  , 


,  .  l-l/2n  -</ 2  ,  l-l/2n. 

y2(r,-ir)  = -ar  n  +o(r  ) 


>  (5.2) 


Thus,  to  dominant  order,  the  two  crack  faces  are  carried  into  the 
curve  represented  by 


^The  shape  of  the  deformation  image  of  the  crack  faces  is  of  course  in¬ 
variant  under  a  rotation  such  as  (3.96). 


i 
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2n  -  1 
2n  +  1 

,  y1<0,  K(n)  *0  .  (5.3)1 

From  (5.1)  it  follows  that  points  near  the  crack-tip  in  the  undeformed 

body  lie  to  the  right  of  this  curve  after  deformation.  We  note  on  the 
2 

basis  of  (5.3)  and  (3.34)  that  the  crack  is  bound  to  open  when  n>3/2, 
regardless  of  the  magnitude  and  nature  of  the  particular  loading  at  in¬ 
finity.  This  conclusion,  which  is  readily  shown  to  hold  also  for 
l<ns3/2,  is  in  marked  contrast  to  the  predictions  of  the  linearized 
theory  for  a  Mode  II  loading. 

It  should  be  emphasized,  however,  that  the  foregoing  result  is  con¬ 
tingent  upon  the  applicability  of  the  Ansatz  (3.4) ,  in  which  the  functions 
are  not  permitted  to  vanish  identically.  For  example,  in  the  presence 
of  a  loading  of  uniform  uniaxial  tension  parallel  to  the  crack  faces,  the 
ensuing  deformation  field  is  homogeneous  and  the  crack  fails  to  open.  In 
this  instance  the  deformation  field  is  non-singular  at  the  crack-tips,  so 
that  (3.4)  cannot  possibly  apply.  We  shall  remark  on  the  question  of  the 
scale  of  the  crack  opening  later  in  this  section. 

Note  that  the  local  deformation  field  (3.34)  may  be  regarded  as  the 
result  of  two  successive  mappings:  a  deformation  of  the  undeformed  body, 
in  which  the  crack  opens,  at  least  when  n>l,  followed  by  a  rigid  rota¬ 
tion,  with  Q  given  by  (3.96).  In  view  of  the  principle  of  objectivity 
ap 

we  shall  therefore  be  entitled  to  base  the  computation  of  the  local  stress 

^The  specific  form  of  K(n) ,  which  is  expressible  in  terms  of  the  Gamma 
function,  is  of  no  particular  interest  here. 

2 

Recall  from  (3.18),  (3.20)  that  the  amplitude  parameter  "a"  cannot  vanish 
in  this  instance. 
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distribution  on  (3.102),  (3.103),  (3.104). 

We  cite  next  the  symmetric  deformation  fields^-  appropriate  to  a 
Mode  I  loading  at  infinity.  Since  U  in  (3.34)  is  an  odd  function,  one 
has  a^=0,  a^  *  a.  This  specialization  leads  to  the  following  results. 
For  l/2<n<  3/2, 


Yl(r,e)  =  bo[U(0)]2r2-1/n  +  o(r2-1/n)  , 

y2(r,0)  =  aU(0)r1  1^2n+o(r2  1^n)  , 

while  for  n>3/2, 

7l(r,0)  *-ir1+1/2nH(0)  +  o(r1+1/2n)  , 

y2(r,0)  =  ar1_1/2nU(0)  +  X(e)r1+1/2n+  o(r1+1/2n) 


I 


(5.4) 


>  (5.5) 


J 


in  which  x  is  furnished  by  (3.90)  with  c2  =  0.  Finally,  for  n=3/2. 


y1(r,0)  =  bQU2(e,3/2)r4//3log  r  +  0(rii//3) 
y2(r,0)  =  ar2/^3U  (0,3/2)  +  o(r2^3) 


•\ 


>(5.6) 


It  is  not  possible,  however,  to  specialize  the  asymptotic  results 
for  the  deformation  so  as  to  arrive  at  an  antisymmetric  local  displace¬ 
ment  field,  at  least  when  n>l.  To  see  this,  recall  that  u  *y  -x  , 

a  J  a  a 

and  observe  that  for  n>l  one  has  m'>l,  whence  (3.34)  give 

^Observe  that  the  associated  displacement  field  ua  =  ya  -  xQ  is  also 
bound  to  be  symmetric  by  virtue  of  the  parity  of  . 
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u^r.0)  -  a1r1"1/2nU(0)  -  r  cos  6  +  o(r) 
u2(r,0)  -  a2r1-1/2nU(0)  -  r  sin  0  +  o(r) 


The  constants  in  (5.7)  evidently  cannot  be  chosen  to  yield  a  displace¬ 
ment  field  of  the  desired  parity. 

As  we  observed  in  the  previous  section,  the  local  asymptotic  solu¬ 
tion  obtained  for  the  Mooney-Rivlin  material  (n «  1)  also  fails  to  admit 
an  antisymmetric  displacement  field.  The  foregoing  conclusions  concern¬ 
ing  nil  lend  support  to  the  conjecture  that  there  does  not  exist  an 
antisymmetric  global  solution  for  all  values  of  the  hardening  parameter. 

At  this  stage  we  calculate  the  associated  local  true-stress  field, 
as  well  as  the  strain-energy  density.  On  account  of  (1.7),  (1.12)  one  has 


t  . *  2W'  (I)F  Fd  -  p6  „ 
a0  ay  0y  a6 


(5.8) 


Suppose  first  that  l/2<n<3/2.  In  this  case,  (3.102)  and  the  first  of 
(3. 120) .together  with  (5.8),  yield 

T11(r,0)*r1  1/,nT11(0,n)  +  o(r1  1^n)  , 

x12(r,0)  *  r_1/,2nT12(0,n)  +  o(r-1/,2n)  ,  ►  (5.9)1 

T22(r,0)  -  r-1T22(0,n) +  o(r-1) 

at 

If  n>3/2,  the  appropriate  stresses  obey 


^Observe  that  the  true  stresses  here  are  referred  to  the  material  polar 
coordinates  (r,0). 
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T11(r,6)  »  o(r~1+1^n)  , 


T12(r,0)  =  r‘1+1/nf12(0,n)  +o(r'1+1/n)  , 


T22(r,0)»r  1T22(0,n)  +  o(r  1) 


>(5.10) 


J 


The  auxiliary  functions  T  appearing  in  (5.9),  (5.10)  are  defined  by 


Tn(0,n)  =  8AbQa^n  2m2nn2  n[u>(8,n)  +  kcos  0]3  cos2£q 
T12(0,n)  =  4Aboa2n'1m2nn3/2_n[a>(0,n)  +  kcos  0]"1/2n  cos  ZQ 


'N 


,>(5.11) 


T22(0,n)  *  2Aa2nm2nn*"  n[w(0,n)  +  iccos  0]  1  ,  J 

with  cos  Cq  given  by  (3.59),  while  T^2  involves  the  hypergeometric 
function  and  will  not  be  cited  explicitly. 

In  the  transition  case  n**3/2,  equations  (4.11),  (4.12),  (4.16) 
may  be  used  to  deduce  the  corresponding  stress  field.  We  note  that  in 
this  instance  t22  is  given  by  the  last  of  (5.9)  or  (5.10)  with  n-3/2 
and  thus  the  representation 


t22^r»9) 


r  lT22^6,n^  +  °<r_1> 


(5.12) 


is  valid  for  l/2<n<a>. 

The  components  of  stress  in  (5.9),  (5.10)  are  those  associated 
with  the  rotated  deformation  fields  (3.102),  (3.103).  If  x*  are  the 
true-stress  components  associated  with  the  original  field  (3.34),  one  has 
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Q  Q  0t 
pa  vB  pv 


(5.13) 


In  which  Q  are  the  components  (3.96)  of  the  rotation  tensor  £. 

Observe  that  a  more  satisfactory  estimate  for  r  ^  than  that  in 
the  first  of  (5.10)  is  precluded  by  the  weak  estimate  for  the  pressure 
in  (3.120).  It  is  interesting  to  recall  in  this  connection  that  the 
asymptotic  solution  in  [3],  [4]  of  the  Mode  I  crack  problem  for  a  class 
of  compressible  elastic  solids  also  fails  to  yield  a  non-degenerate 
estimate  for  the  stress  at  sufficiently  large  values  of  the  cor¬ 

responding  hardening  parameter.  The  most  singular  of  the  stress  com¬ 
ponents  in  (5.9),  (5.10)  is  T22’  becomes  infinite  at  the  crack-tip 

like  r  According  to  (5.11),  T22>0,  an{*  hence  T22  is  tensile  ^or 

all  small  enough  values  of  r. 

We  now  determine  the  dominant  character  of  the  true  stresses  when 
the  latter  are  referred  to  the  spatial  coordinates  y  .  Since  such  a 
representation  of  the  stresses  depends  on  the  availability  of  an  in¬ 
vertible  estimate  for  the  local  deformation  we  confine  our  attention  at 
present  to  nil.  One  draws  from  (3.34)  that 


(2n-l)/2n  .  ,  -  -  . 

r^.  ~r  if  r^.  *  y1  +  y 


.2. 

L"*  ' 


(5.14) 


Equations  (5.9),  (5.14)  imply  that  for  l^n<3/2, 


2(n-l)N 

2n-l 


lll 


(5.15) 


For  n>3/2,  (5.10)  yield 
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-2(n-l)N 

2n-l 


-2(n-l)' 

2n-l 


-2n 

2n-l 


11 


12 


22 


(5.16) 


We  observe  that  the  singularity  of  t^2  in  (5.15),  (5.16)  becomes  pro¬ 
gressively  more  severe  with  increasing  values  of  the  hardening  parameter 
and  —  for  the  range  of  n  under  consideration  —  is  always  stronger  than 
that  predicted  by  the  linear  theory. 

From  (1.18),  (3.34)  we  gather  that 


W(I)  »Aa2nGn(e)r-1+o(r-1),  r-0,  n>l/2  ,  (5.17) 


so  that  the  strain-energy  density  has  a  singularity  near  the  crack-tips 
which  is  of  the  same  order  as  its  counterpart  in  the  linear  theory. 

Finally  we  supply  an  estimate  for  the  amplitude  parameter  "a", 
valid  for  small  loads  at  infinity.  In  this  connection  we  recall  a  conser¬ 
vation  law  appropriate  to  the  equilibrium  theory  of  finite  plane  strain: 
for  every  curve  C  that  is  the  boundary  of  a  finite  regular  subregion 
of  .n , 

(Wn  -  FTs)dfc  =  0  ,  (5.18)1 

/w  rw  rv 

C 

provided  W  and  F  are  the  strain-energy  density  and  the  deformation 
gradient  tensor,  whereas  £  is  the  unit  outward  normal  of  C  and  s 
is  the  Piola  traction, 


(5.18)  is  an  adaptation  to  plane  strain  of  a  conservation  law  originally 
due  to  Eshelby  [16].  The  importance  of  this  conservation  law  in  fracture 
mechanics  was  first  recognized  by  Rice  [17]. 
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s  ■  o  n 


on 


(5.19) 


The  conservation  law  (5.18)  holds  for  both  compressible  and  Incompressible 
materials,  and  is  readily  confirmed  with  the  aid  of  the  divergence  theorem 
and  the  appropriate  field  equations. 

For  the  present  purpose  we  consider  first  the  Mode  I  loading  case 
and  assume  that  the  solution  to  the  nonlinear  problem  for  sufficiently 
small  loads  is  uniformly  approximated  by  the  corresponding  solution  from 
the  linear  theory,  at  all  material  points  a  finite  distance  from  the  crack- 
tips.  A  limit  process  depending  on  this  assumption  enables  one  to  deduce  from 
(5.18)  and  from  the  available  solution  of  the  linearized  problem  the  de¬ 
sired  small-load  amplitude  estimate.  Since  the  required  calculation  is 
strictly  parallel  to  one  spelled  out  in  detail  in  [3],  we  cite  the  en¬ 
suing  estimate  directly: 


o  2.  n-2  „ 

-2n-  a  bn  +o(o2) 


4pAm 


2n-l 


as  a-*  0 


(5.20) 


Here  a  is  the  applied  normal  stress  at  infinity,  2b  is  the  crack 
length,  while  p  is  the  shear  modulus  at  infinitesimal  deformations. 

According  to  a  conclusion  reached  earlier  in  this  section,  the 
finite  theory  —  in  contrast  to  linearized  elastostatics  —  predicts  that 
the  crack  opens  in  the  vicinity  of  its  tips  also  for  a  Mode  II  loading, 
at  least  for  a  certain  range  of  the  hardening  parameter.  Clearly,  if  the 
solution  to  the  linearized  Mode  II  problem  is  to  approximate  its  counter¬ 
part  in  the  finite  theory  uniformly  at  small  loads  on  every  bounded  and 
closed  material  point  set  that  excludes  the  crack-tips,  the  lengths  of 
the  opening  crack-segments  would  have  to  tend  to  zero  with  the  load 
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intensity.  Making  this  assumption,  invoking  the  conservation  law  (5.18) 
once  again,  and  appealing  to  the  Mode  II  results  of  the  linear  theory 
summarized  in  (2.12),  (2.13),  (2.14),  one  is  led  to  precisely  the  esti¬ 
mate  obtained  for  Mode  I: 


,  2.  n-2  . 

.  _  2n  o  bn  ,  /  2. 

L  =  a  = - —  r  +  o(o  ) 


4yAm 


2n-l 


(5.21) 


in  which  a  is  now  the  shear  stress  applied  at  infinity.  The  length 
parameter  L,  in  view  of  (5.3),  governs  the  scale  of  the  crack-opening, 
which  is  thus  found  to  be  a  second-order  effect  at  small  shearing  loads. 
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Appendix:  Representation  of  functions  biharmonic  in  a  neighborhood 
of  infinity. 

In  this  paper  a  function  U  is  called  biharmonic  on  an  open 

plane  region  ft,  which  need  not  be  simply  connected,  if  it  is  four  times 

4 

continuously  differentiable  on  ft  and  V  U  =  0  on  ft. 

Let  71  be  a  neighborhood  of  infinity  in  the  (x^.x^) -plane,  so 
that  7\  is  the  exterior  of  some  circle  centered  at  the  origin,  and  let 
7 be  the  simply-connected  domain  obtained  by  deleting  from  7(  its 
intersection  with  the  positive  x^-axis.  As  shown  by  Muskhelishvili 
[14]  in  his  analysis  of  the  structure  of  the  Airy  stress  function  on 
multiply  connected  domains,  any  function  U  biharmonic  on  7('  ,  whose 
partial  derivatives  of  the  second  and  higher  orders  admit  continuous  ex¬ 
tensions  to  7 <,  may  be  represented  in  the  form 


U(x1,x2)  -  Re  [z'f(z)  +  x(z)}  , 


(1) 


where 


^(z)  -  Azlog  z  +  Blog  z  +  4'*(z)  , 

X(z)  =  Czlog  z  +  Dlog  z  +  x*(z) 


^  (2) 


Here  z  ■  x^ +  ix2  ■  pe^  ,  log z  refers  to  the  principal  branch  of  the 
logarithm,  while  V*(z)  and  x*(z)  are  functions  analytic  on  % 

Let  Ap  B^,  C^,  and  A2,  B2>  C2,  D2  denote  the  real  and 
imaginary  parts  of  A,  B,  C,  D,  respectively.  Invoking  the  Laurent  ex¬ 
pansions  on  71  of  Y*(z)  and  x*(z),  one  readily  deduces  from  (1)  and 


i  ■ 


(2)  the  representation 


2  2 

U(x1,x2)-A1p  logp-A2p  9 

+  [  (B^  +  C^)cos  9+  (B2  +  C2)sin  9]plog  p 
+  [  (B^  -  C2)  sin  9  -  (B2  +  C-^cos  9]p9 

GO 

+  D^logp-D29+  £  [a^cos  j9  +  bjSin  j9 

j“  -  CD 

+  CjCos(j  -  2)9+ djSin(j  -  2>9]p  ,  (3) 

valid  on  7('  ,  the  coefficients  a^ ,  b  ^  ,  c^ ,  d^  being  real  constants. 

Suppose  now  U  is  biharmonic  on  71  and  hence  .continuous  on  71 
together  with  all  its  partial  derivatives  up  to  the  fourth  order.  Then 
(3)  demands  that 

A2=0;  B1-C2»0;  B^CjjO,  D2  -  0  .  (A) 

Suppose  further  that 

U(x^,x2)  ■  0(1)  as  p-*oo  .  (5) 

This  requirement  evidently  yields  the  additional  restrictions 
»j-0.  Bj  +  Cj-O,  B2  +  C2-°,  Dj-O,  VVW°  for  j<0  .(6) 

In  view  of  (A),  (6),  the  representation  (3)  finally  reduces  to 


A 
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U(x^,x2)=^  [a^cos  jcp  +  bjSin  jcp+ CjCos(j  -  2)cp+ d^sin(j  -  2)cp]p  ^  .  (7) 

j-o 

This  expression  is  valid  on  7\ ,  the  infinite  series  being  uniformly  and 
absolutely  convergent  on  every  closed  subset  of  71 .  Thus  any  function 
biharmonic  on  7[  and  bounded  at  infinity  necessarily  admits  the  re¬ 
presentation  (7)  . 

Acknowledgment 

The  author  wishes  to  express  his  appreciation  to  Professor 
Eli  Sternberg  for  his  invaluable  guidance  throughout  this  investigation. 
I  am  also  greatly  indebted  to  Professor  J.K. Knowles  for  most  helpful 


discussions. 


-63- 


References 


[1]  Rice,  J.R.,  Mathematical  analysis  in  the  mechanics  of  fracture. 
Fracture,  Vol.II,  Chapter  3,  Academic  Press,  New  York,  1968. 

[2]  Wong,  F.S.  and  Shield,  R.T.,  Large  plane  deformations  of  thin 
elastic  sheets  of  neo-Hookean  material,  Zeltschrift  ftlr  angewandte 
Mathematik  und  Physik,  20  (1969)  2,  p.176. 

[3]  Knowles,  J.K.  and  Sternberg,  E.,  An  asymptotic  finite-deformation 
analysis  of  the  elastostatic  field  near  the  tip  of  a  crack.  Journal 
of  Elasticity,  3  (1973)  2,  p.67. 

[4]  Knowles,  J.K.  and  Sternberg,  E. ,  Finite-deformation  analysis  of  the 
elastostatic  field  near  the  tip  of  a  crack:  Reconsideration  and 
higher  order  results.  Journal  of  Elasticity,  4^  (1974)  3,  p.201. 

[5]  Knowles,  J.K.,  The  finite  anti-plane  shear  field  near  the  tip  of 
a  crack  for  a  class  of  incompressible  materials.  International 
Journal  of  Fracture,  13  (1977)  5,  p.61. 

[6]  Lo,  K.K.,  Finite-deformation  crack  in  an  infinite  body  under  anti- 
plane  simple  shear ,  International  Journal  of  Solids  and  Structures, 

13  (1977),  p.1219. 

[7]  Knowles,  J.K.  and  Sternberg,  E.,  Discontinuous  deformation  gradients 
near  the  tip  of  a  crack  in  finite  anti-plane  shear:  an  example, 
Journal  of  Elasticity,  H)  (1980)  1,  p.81. 

[8]  Knowles,  J.K.  and  Sternberg,  E.,  Anti-plane  shear  fields  with  dis¬ 
continuous  deformation  gradients  near  the  tip  of  a  crack  in  finite 
elastostatics ,  Technical  Report  No. 43,  ONR  Contract  N00014-75-C-0196, 
September  1979.  To  appear  in  Journal  of  Elasticity. 

[ 9 ]  Knowles ,  J . K . ,  Finite  elastostatic  fields  with  unbounded  deformation 
gradients ,  in  Finite  Elasticity,  Applied  Mechanics  Division,  ASME, 
Vol.27  (1977),  p. 23. 

[10]  Knowles,  J.K.,  Crack  problems  in  finite  elastostatics.  Proceedings, 
Symposium  on  Mathematical  Problems  in  Fracture,  AMS  and  SIAM, 

New  York,  March  1978. 

[11]  Knowles,  J.K.,  On  some  inherently  nonlinear  singular  problems  in 
finite  elastostatics.  Proceedings,  Eighth  U.S. National  Congress  of 
Applied  Mechanics,  UCLA,  June  1978. 


-64- 


[12]  Sternberg,  E.,  Some  recent  advances  In  the  application  of  nonlinear 
elastostatlcs  to  singular  problems,  in  Trends  in  Solid  Mechanics 
1979,  Proceedings  of  the  Symposium  dedicated  to  the  65th  Birthday 
of  W.T.Koiter,  Delft  University  Press,  Sijthoff  and  Noordhoff 
International  Publishers,  1979. 

[13]  Abeyaratne,  R.C.,  Discontinuous  deformation  gradients  in  plane 
finite  elastostatlcs  of  incompressible  materials.  Technical  Report 
No. 41,  ONR  Contract  N00014-75-C-0196 ,  June  1979.  To  appear  in 
Journal  of  Elasticity. 

[14]  Muskhelishvili,  N.I.,  Some  basic  problems  of  the  mathematical  theory 
of  elasticity.  Second  English  edition,  P. Noordhoff  Ltd.,  Groningen, 
1963. 

[15]  Abromowitz,  M.  and  Stegun,  I. A.,  eds..  Handbook  of  mathematical 
functions,  National  Bureau  of  Standards,  Applied  Mathematics  Series 
55,  1964. 

[16]  Eshelby,  J.D.,  The  continuum  theory  of  lattice  defects.  Solid  State 
Physics,  Vol.3,  Academic  Press,  New  York,  1956. 

[17]  Rice,  J.R.,  A  path-independent  integral  and  the  approximate  analysis 
of  stress  concentration  by  notches  and  cracks.  Journal  of  Applied 
Mechanics,  _35  (1968)  2,  p.379. 


[«1«] 


x2 


FIGURE  2.  SLAB  WITH  CRACK  AND  COORDINATES. 
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